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Prediction  of  transverse-mode  selection  in  double  heterojunction  lasers  by 
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Jerome  K.  Butler 
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(Received  6  July  1982;  accepted  for  publication  23  September  1982) 

Transverse-mode  selection  is  characterized  for  GaAs/AlGaAs  double  heterojunction  lasers  from 
optical  field  and  electron/hole  interaction.  The  electron/hole  distribution  determined  from  a 
solution  of  the  ambipolar  diffusion  equation  provides  the  necessary  information  about  gain/mode 
coupling  to  predict  the  current  at  threshold.  Lasing  power  out  versus  current  solutions  provide 
information  about  internal  differential  quantum  efficiency.  Theory  is  matched  to  experiment  for  a 
multimode  laser  with  one  heterojunction  having  a  very  small  index  step.  It  is  found  that  the  laser’s 
characteristics  over  a  temperature  and  current  range  are  predicted  by  adjusting  the  active-layer 
refractive  index  as  determined  from  far-field  measurements. 


PACS  numbers:  42.55.Px,  42.60.By 

I.  INTRODUCTION 

In  recent  years,  there  has  been  considerable  effort  spent 
fabricating  various  semiconductor  laser-device  geometries 
for  control  of  mode  operation.  In  GaAs/AlGaAs  devices, 
the  transverse-mode  operation  is  usually  governed  by  the 
growth  of  various  layer  thicknesses  and  index  steps  at  the 
grown  heterojunctions.  Double  heterojunction  lasers  which 
confine  the  optical  field  to  the  thin  active  layers  usually  oper¬ 
ate  in  the  fundamental  mode.  Furthermore,  the  lasing  mode 
is  independent  of  drive  because  small  index  changes  due  to 
injection  have  negligible  effects  on  mode  selection.  While 
transverse-mode  lasing  can  be  well  controlled  by  epitaxial 
heterojunctions,  lateral-mode  operation  in  contemporary  la¬ 
sers  is  drive  sensitive.  Stripe-contact  devices1-3  have  lateral 
guides  defined  only  by  the  current  density  which  affects  the 
gain  distribution  while  more  sophisticated  structures4-7 
have  grown  layers  whose  thicknesses  are  functions  of  lateral 
positions.  Thus,  in  these  devices  the  lateral  modes  are  shaped 
by  the  lateral  effective  index  as  well  as  gain  variations. 

Because  the  gain  distribution  plays  a  major  role  in  de¬ 
fining  mode  stability  in  lasers,  it  is  important  to  understand 
the  mechanisms  of  electron/hole  transport  in  active  layers. 
In  this  paper,  we  discuss  the  transverse-mode  operation  of 
lasers  which  are  affected  by  twe  electron/hole  distribution. 
For  the  first  time  the  ambipolar  diffusion  equation  is  solved 
for  the  active  layer  in  an  injection  laser.  The  direct  conse¬ 
quence  of  this  solution  is  the  ability  to  predict  crossover  of 
competing  modes  in  the  active  layer.  In  contrast  to  ambipo¬ 
lar  diffusion,  assuming  simple  electron  injection  does  not 
explain  this  phenomena  in  sufficient  detail. 

In  this  work  we  investigate  mode  selection  from  the 
standpoint  of  the  exact  selection  process.  This  includes  mode 
gain  coupling  and  differential  quantum  efficiency.  We  con¬ 
centrate  on  the  transverse  (perpendicular  to  the  junction 
plane)  modes  in  a  double  heterojunction  laser.  The  excess 
carrier  distribution  is  assumed  ambipolar  and  a  complete 
solution  is  obtained.  From  mode/gain  coupling,  it  is  possible 


to  predict  the  threshold,  and  from  the  slope  efficiency,  it  is 
possible  to  predict  mode  crossover. 

The  laser  we  apply  this  to  is  a  broad-area  double  hetero¬ 
junction  GaAs  device  with  a  small  index  step  at  one  of  the 
heterojunctions.  The  laser  is  heated  over  a  temperature 
range  and  the  first  two  modes  are  investigated.  It  is  found 
that  cavity-mode  preference  is  very  sensitive  to  temperature 
and  current  level.  Application  of  the  diffusion  equation  is 
used  to  describe  the  crossover  of  the  first  two  modes  and  it  is 
established  that  the  mode  with  the  highest  slope  efficiency 
will  dominate.  Far-field  patterns  are  described  by  standard 
multilayer  waveguide  techniques.  Finally,  a  calculation  is 
made  to  estimate  the  effect  of  free  carriers  on  the  refractive 
index  in  the  material. 

II.  THEORY 

A.  Approach 

Excess  carrier  movement  in  the  active  layer  of  an  injec¬ 
tion  laser  is  described  using  the  ambipolar  diffusion  equa¬ 
tion.  Quasi-Fermi  level  continuity  as  well  as  electron  and 
hole  current  continuity  is  imposed  at  each  interface.  The 
quasi-Fermi  level  location  is  calculated  using  a  nonparabolic 
riC  extremum,  and  parabolic  X]C  and  T ,,,,  extrema.  The 
boundary  condition  is  set  up  in  such  a  way  as  to  eliminate  the 
need  for  a  solution  within  the  space  charge  layer  surround¬ 
ing  each  interface.  The  spontaneous  recombination  lifetime 
is  calculated  without  imposing  quasimomentum  conserva¬ 
tion  for  parabolic  bands.  Above  threshold,  a  stimulated  re¬ 
combination  term  is  included  in  the  ambipolar  diffusion 
equation  and  the  gain  is  assumed  to  be  pinned  at  the  thresh¬ 
old  value.  The  lasing  modal  shape  for  the  transverse  modes 
(perpendicular  to  the  metallurgical  junctions)  is  calculated 
by  solving  Maxwell’s  equations  for  a  multilayer  dielectric 
wave  guide. 

B.  Ths  ixcmi  carrier  distribution 

The  diffusion  equation  evolves  along  the  standard 
route.  We  demand,  as  is  customary,  that  each  carrier  type 
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satisfies  the  continuity  equation8 

DpV2p-/ipV-[pE)+gp  -p/rp  =  dp/dt,  (la) 

Dn V2n  +  p„ V-[nE )  +  g„  -  n/r„  =  dn/dt,  (lb) 

where  n,  p  —  total  electron,  hole  carrier  density;  D„,  Dp 
=  electron,  hole  diffusion  constant;  g„,gp  =  electron,  hole 
generation  rate;  and  t„,tp  =  electron,  hole  mean  recombin¬ 
ation  time.  We  further  demand  that  Poisson’s  equation  be 
satisfied: 

V-£  =  (Airp/e),  (2) 

where  p  is  the  space  charge  density  and  e  is  the  dielectric 
constant. 

Rather  than  solve  Eqs.  (1)  and  (2)  exactly,  we  take  an 
approximate  route  where  we  assume  that  the  internal  field 
between  the  electron  and  hole  charge  densities  is  strong 
enough  to  guarantee  charge  neutrality.  Although  this  condi¬ 
tion  is  violated  in  the  space  charge  layer  surrounding  the 
heterojunctions,  we  can  formulate  the  problem  in  such  a  way 
that  only  diffusion  processes  outside  the  heterojunctions 
need  be  considered.  Under  this  assumption,  the  excess  elec¬ 
tron  and  hole  densities,  8n  and  8p,  respectively,  must  be 
equal: 

8n  =  n  —  n0  =p  —p0  =  8p.  (3) 

Here,  n0,p0  designate  the  thermal-equilibrium  carrier  densi¬ 
ties.  If  we  take  {dn0)/(dx)  —  (dp0)/(dx)  =  0,  we  can  substitute 
Eq.  (3)  into  Eqjl)  to  obtain  two  equations  for  the  two  unk¬ 
nowns  8p  and  E.  The  terms  involving  V  E  can  be  eliminated 
by  multiplying  Eq.  ( la)  by  np„  and  Eq.  ( lb)  by ppp  and  add¬ 
ing  to  obtain9 

D *V2[8p)  - p*E-V(8p)  +g'-^  =  -3&L ,  (4) 

r  dt 

where 


(~2x  1018  cm'). 

The  electric  field  in  Eq.  (4)  refers  to  both  the  internal 
and  applied  fields.  Under  the  conditions  encountered  in  an 
injection  laser,  these  fields  exert  only  a  small  direct  influence 
on  carrier  distribution,  and  as  such  can  be  ignored  in  Eq.  (4). 
In  other  words,  the  electric  field  does  not  alter  the  actual 
excess  carrier  spatial  distribution,  but  rather  only  causes  car¬ 
riers  to  drift.  This  is  the  case  if  the  drift  length  is  much  less 
than  the  diffusion  length,  i.e., 

P*Et<{D*t)'12.  (5) 

This  condition  is  met  in  the  solutions  ultimately  obtained. 
The  major  effect  of  the  internal  field,  that  of  charge  neutra¬ 
lity,  is  still  accounted. 

If  we  assume  the  only  excess  carrier  drain,  other  than 
spontaneous  recombination  8p/r,  is  stimulated  recombina¬ 
tion,  then  for  an  infinitely  long  cavity  and  time-averaged 
values10 

r.-S&L.llB*).  (6, 

dt  rjs  hv 

where  P{x)  =  time-averaged  power  density,  g  =  gain,  i js 
=  stimulated  quantum  efficiency,  and  hv  =  photon  energy. 
For  100%  carrier  confinement,  t]%  is  the  internal  differential 
quantum  efficiency.  In  these  calculations  we  put  tjs  =  1. 
Equation  (4)  becomes 

D  *V2(<5p)  —  (8p)/r  =  r,  (7) 

where/)  *  =  ambipolar  diffusion  constant,  r  =  spontaneous 
recombination  lifetime,  and  r  =  stimulated  recombination 
rate.  The  value  r  is  equated  to  zero  in  the  region  where  g  <  0. 
The  spontaneous  recombination  lifetime  is 

T=l/lB[n0+p0  +  8p)],  (8) 

where  B  is  1.3  X  10~ 10  cmVsec. 11,12 


8p  —  excess  carrier  density, 
np„D„+pppDp 


Wn  +  PPp 
np„  +PPP 


*  =  p-jb^o-PA 


g  =g'n  =  gp, 


8p  _Po  +  dp  Pq  _  n0  4-  8n  n0 

T  tp  tp„  r-  r»„  ’ 

g'*<8p  —  excess  carrier  generation  rates,  and 

r  =  excess  carrier  recombination  lifetime. 


We  have  implicitly  assumed  that  the  net  recombination  rate 
(8p)/r  is  the  same  for  electrons  as  for  holes,  as  it  must  be.  We 
have  further  assumed  that  the  equilibrium  generation  rate  is 
equal  to  the  equilibrium  recombination  rate  for  both  elec¬ 
trons  and  holes. 

For  the  ambipolar  diffusion  Eq.  (4),  observe  that  the 
ambipolar  diffusion  rate  for  excess  carriers  is  determined  by 
the  average  value  of  D„  and  Dp ,  weighted  by  np„  and  ppp .  In 
other  words,  the  two  oppositely  charged  clouds  interact  with 
each  other  through  the  internal  field  in  such  a  way  as  to 
satisfy  Poisson’s  equation  and  must  then  diffuse  at  a  natural¬ 
ly  acceptable  rate.  This  is  crucial  because  of  the  very  large 
injected  charge  densities  encountered  in  injection  lasers 


C.  The  boundary  conditions 

At  each  boundary  within  the  device  there  exists  a  space 
charge  layer  with  a  width  severely  reduced  under  conditions 
of  forward  bias.  Because  quasicharge  neutrality  is  violated 
within  this  layer,  we  match  across  the  layer  under  the  as¬ 
sumptions  it  is  is  thin  enough  and  internal  fields  are  strong 
enough  that  carriers  are  swept  through  without  significant 
recombination.  Hence,  we  demand  that  the  electron  and 
hole  current  densities  be  the  same  on  each  side  of  the  space 
charge  layer.  These  densities  are  written 

Jn  =  qnp„E  +  qD„  V[8p),  (9a) 

lp  =  qpppE  -  qDpV{8p).  (9b) 

Here,  £_is  the  total  field  and  q  is  the  electron  charge.  Al¬ 
though  E  can  be  ignored  in  Eq.  (7),  where  it  has  a  small  effect, 
it  must  be  included  in  the  boundary  conditions  because  both 
drift  and  diffusion  terms  are  significant. 

In  addition  to  continuity  of  current  densities,  we  also 
impose  continuity  of  hole  and  electron  quasi-Fermi  levels 
across  each  layer  interface.  Although  it  is  clear  this  align¬ 
ment  of  Fermi  levels  implies  band  bending,  i.e.,  space 
charge,  we  assume  an  abrupt  heterojunction  and  ignore  such 
effects. 
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D.  The  electromagnetic  field  distribution 

The  standard  multilayer  slab  dielectric  waveguide 
model  is  used  to  describe  the  transverse  far-field  pat¬ 
terns.  1  * 14  Standard  calculations  yield  propagation  constants 
Y  =  a  +  j[3  and  fields.  For  instance,  for  a  TE  mode,  we  can 
express  the  rwth-mode  field  distribution  for  the  ;th  layer  as 

E |.t|  =  E,n  cos (h  '„x  +  6  I,  (10) 

where  h  is  the  complex  eigenvalue  and  x  is  the  coordinate 
across  the  active  region.  The  overall  field  solution  is  obtained 
by  matching  the  field  components  at  each  interface.  The  far- 
field  pattern  then  becomes  proportional  to  the  Fourier  trans¬ 
form  of  the  near  field  on  a  lasing  facet. 


(L*)2 


i=  1,3, 


E.  The  diffusion-equation  solution 

Assume  a  three-layer  device  with  the  central  layer  ac¬ 
tive  (see  Fig.  3).  Also,  assume  the  mode  is  propagating  down 
the  z  axis  with  x  perpendicular  to  the  guide.  The  excess  car¬ 
rier  diffusion  equations  for  this  case  become 

d  2(<Sp) 
dx- 

d2(6p)  Sp  g(x)PJx) 
dx2  (L  *)2  VM 
where  L*  =  \D  *V)l/\ g{x)  is  the  localized  gain  and  P„  (x)  is 
the  optical  Poynting  vector  component  along  the  propaga¬ 
tion  direction. 

These  equations  are  linearized  by  replacing  g(x)  with  its 
weighted  average  within  the  active  region10 

f-  ,*(*)£ U*)dx  gmNm 


i  —  2, 


(Ha) 


(Hb) 


g  = 


(12) 


SE2m(x)dx  Tm 

where  Nm  is  the  field  normalization  constant  and  Tm  is  the 
mode  intensity  integrated  over  the  active  layer.  It  should  be 
noted  that  g(x)  =  0  outside  the  active  region.  Further,  as¬ 
suming  no  losses  in  layers  external  to  the  active  region,  gm 
represents  the  propagation  gain  of  the  mth  mode.  With  D  * 
and  r  assumed  to  be  piecewise  constant  layer-to-layer,  Eq. 
(11b)  becomes 


d2{Sp) 
dx 2 


A, 


(£?)2 


cos2(&^x-M; 


(13) 


where 


A- 


-J. 


D*Tmhvr,, 
PAx\dx. 


Ignoring  the  minor  effect  of  the  layer  3-4  interface  on  the 
carrier  distribution,  the  solution  to  the  diffusion  equation  in 
the  various  layers  becomes 

Sp3[x)  =  <?3  cosh {x/L  ?)  +  F3  sinh (x/L  ?),  (14a) 

Sp2  =  G2  cosh(x/I  J)  4-  F2  sinh(x/I  ?) 

+  («*) 


&PM)  =  Ci  cosh(x/L  f)  +  F,  sinh(x/L  f ). 


(14c) 


With  these  eigenfunctions,  the  solution  follows  directly 
from  matching  the  current  density  and  quasi-Fermi  levels  at 
each  interface.  The  Fermi-level  continuity  condition  is  ap¬ 
plied  using  Kane’s  nonparabolicity  for  the  f1(  minimum 
and  a  parabolic  expression  for  both  r]3V  and  the  Xlc  mini¬ 
ma.  The  effect  of  aluminum  concentration  on  carrier  mobil¬ 
ity,  effective  mass,  and  band  gap  is  included.  The  volume 
recombination  rate  is  related  to  the  optical  gain  of  the  laser  in 
linear  fashion  using  coefficients  obtained  by  integrating  over 
parabolic  bands  without  quasimomentum  conservation.  The 
recombination  lifetime  is  found  in  a  similar  manner.  The 
diffusion  coefficients  are  calculated  using  the  generalized 
Einstein  equation.  The  complete  formulation  of  the  problem 
is  presented  elsewhere.15  The  necessary  roots  are  calculated 
on  a  computer  using  a  real-root  searching  routine.  The  final 
solution  provides  current-voltage  curves,  power  output-cur- 
rent  curves,  threshold  current  density  of  each  mode,  and 
finally  electron  and  hole  injection  efficiencies  into  the  active 
'layer. 

F.  Validation  of  numerical  results 

This  theory  has  been  applied  to  double  heterojunction 
lasers,  single  heterojunction  lasers,  and  separate  electron- 
optical  confinement  lasers.  Close  agreement  has  been  found 
with  experimental  values.  Two  typical  results  for  a  double 
heterojunction  laser  are  shown  in  Figs.  1  and  2.  The  first 
shows  threshold  current  density  as  the  active  region  thick¬ 
ness  is  reduced,  and  the  second  shows  threshold  as  a  function 
of  the  aluminum  concentration  in  the  passive  regions.  The 
slope  of  Fig.  1  is  5.2  kA /cm2 /im,  a  value  which  compares 
closely  with  available  data.1617 


>- 


FIG.  1.  Calculated  hole  injection  efficiency,  relative  peak  power,  and 
threshold  current  density  of  a  symmetric  DH  laser  as  a  function  of  cavity 
width.  The  active  layer  doping  isp- 1 0 1 7  cm  ~ 3  and  the  passive  a -layer  doping 
is  10'*  cm-3.  The  passivep  layer  hasaSx  10"  cm" 3  doping  level.  The  AlAs 
percentage  in  the  outer  layers  is  20%. 
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FIG.  2.  Calculated  threshold  current  density  as  a  function  of  AlAs  percen¬ 
tage  in  the  GaAlAs  passive  p  layer  of  a  symmetric  DH  laser  with  a  l-ytm 
wide  GaAs  (p-1017  cm-3)  active  layer.  Three  passive p-layer  doping  levels 
are  shown.  The  passive  n-layer  doping  level  is  1 X  1018  cm-3. 


III.  EXPERIMENT 

The  lasers  characterized  in  this  experiment  are  broad- 
area  double  heterojunction  lasers  grown  by  liquid  phase  epi¬ 
taxy  methods.  Figure  3  gives  the  geometry  of  one  such  laser, 
PL174-1-31.  Other  lasers  of  this  type  were  characterized, 
but  PL  174- 1-31  was  selected  for  its  transverse  modal  behav¬ 
ior  with  temperature  due  to  the  small  refractive  index  step  at 
one  heterojunction.  Estimating  the  refractive  index  of  the 
aluminum  layers  is  accomplished  using  experimental 
data18-19  at  a  lasing  wavelength  A  =  9000  A.  This  gives 
At)  =  0.62x,  where  x  is  the  fraction  of  aluminum  in  the  solid. 
The  laser  has  a  cross  section  of  214x356  /im,  where  the 
latter  is  the  cavity  length.  The  chip  thickness  is  77  fxm.  The 
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FIG.  3.  Geometry  of  PL174-I-31.  For  the  excess  carrier  calculation,  only 
layers  1, 2,  and  3  are  used. 


laser  is  mounted p-side  up  on  the  header  and  has  no  reflective 
coating  on  the  facets. 

The  device  is  driven  with  250-nsec  pulses  with  a  0.01  % 
duty  cycle  to  minimize  internal  heating.  While  driven,  the 
laser  is  heated  over  temperature  range  of 20-40  °C.  The  tem- 
perature/current-drive  characteristics  of  the  laser  are  char¬ 
acterized  by  its  far-  and  near-field  dependence. 

To  measure  the  far  field,  the  laser  is  placed  on  a  rotating 
platter  with  an  axis  turned  by  a  2-rph  timing  motor.  The 
laser's  output  is  collected  by  monochromator  and  amplified 
by  an  RCA  7102  photomultiplier  tube  (PMT).  Typical  mea¬ 
surement  parameters  are  2-A  slit  width  and  750-V  bias  for 
the  PMT.  The  output  of  the  PMT  is  fed  to  a  lock-in  amplifier 
triggered  by  a  reference  signal  from  the  pulser.  An  HP  7560 
A  log  converter  and  HP  2470  AX  10  dc  multiplier  process 
the  signal  for  they  axis  of  a  chart  recorder.  An  impedance 
matching  network  allows  they  axis  to  be  calibrated  to  three 
decibels  per  inch  deflection.  The  calibration  is  good  for  over 
five  decades  of  response.  The  speed  of  the  x  axis,  combined 
with  the  2-rph  timing  motor  presents  24°  rotation  per  inch 
deflection  for  the  x  axis.  An  alternative  presentation  is  to 
load  the  PMT  with  a  50-/2  terminator  and  display  the  output 
pulse  on  a  sampling  oscilloscope  with  a  50-/2  input.  This 
provides  a  time  decomposition  of  the  light-output  pulse  for 
ready  comparison  to  the  current-input  pulse.  Patterns  were 
typically  made  with  J  ~  1.  l/th . 

The  measurement  scheme  for  the  near  field  is  accom¬ 
plished  by  attaching  the  laser  mount  to  a  micropositioner 
which  has  a  three-dimensional  adjustment.  A  microscope 
objective  with  160  magnification  is  mounted  opposite  the 
laser  facet  and  a  video  camera  collects  the  output  of  the  ob¬ 
jective  lens  and  displays  the  image  on  a  monitor. 

Peak  threshold  current  values  for  laser  PL  174- 1-3 1  are 
shown  in  Fig.  4.  Triangles  denote  the  fundamental  mode 
whereas  crosses  mark  the  second.  In  the  heat-sink  tempera¬ 
ture  range  20-28  °C,  only  the  second  mode  is  present;  for  28- 
42  °C,  both  fundamental  and  second  modes  are  present;  and 
above  42  °C,  only  the  fundamental  mode  propagates.  The 
lasing  wavelengths  for  the  two  modes  are  shown  in  Fig.  5. 

The  multilayer  waveguide  model  is  used  to  describe  the 
transverse  far-field  patterns  of  the  modes.  Figure  6  is  an  ex¬ 
ample  of  a  match  of  the  second  mode  at  20  °C.  The  solid 
curve  is  the  measured  pattern,  and  the  dashed  curve  is  the 
match.  The  bump  on  one  shoulder  at  21.6°  of  Fig.  6  is  the 
result  of  light  coupling  to  the  lossy  substrate.  This  bump  can 
be  accurately  positioned  to  match  the  experimental  pattern 
by  applying  Snell’s  law  at  the  substrate-air  interface.  The 
internal  angle  from  the  normal  to  the  facet  6 '  is  given  by 
tan  0  '~A  4//?~»/4  sin  6.  Here,  h  4  is  the  imaginary  part 

of  h4,  the  real  propagation  constant 0  =  Im{  y|,  6  is  the  ex¬ 
ternal  angle  to  the  normal  to  the  facet,  and  tj4  is  the  refractive 
index  of  the  substrate.  The  relative  intensity  or  amplitude  of 
the  bump  is  adjusted  by  the  separating  width  d}. 

6PP  is  defined  as  the  angle  separating  the  two  major 
lobes  of  the  second  mode.  The  solid  line  in  Fig.  7  shows  how 
8PP  changes  over  the  temperature  range.  The  dashed  line  in 
Fig.  7  is  the  multilayer  6PP  versus  the  active-region  refractive 
index  i/2.  Over  the  entire  range  20-42  “C  the  active-region 
index  depression  necessary  to  match  experimental  evidence 
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FIG.  4.  Peak  threshold  current  vs  heat-sink  temperature.  The  triangles  are 
data  points  Tor  the  fundamental  mode  and  the  crosses  for  the  second  mode 
of  PL174-1-31. 


ANGLE  FROM  NORMAL  (DEG  ) 

FIG.  6.  Experimental  (solid  line)  and  theoretical  (dashed  line)  far-field  pat¬ 
terns  for  the  second  mode  of  PL  174- 1-31  at  20*. 


is  Si)  =  —  0.002.  The  fundamental  mode  also  is  marked  by 
this  refractive  index  decrease  over  the  temperature  range. 
The  half-power  full-width  of  the  fundamental  transverse 
mode  ffl  decreased  over  28-42  °C  in  a  manner  consistent 
with  the  Sij  calculated  for  the  second  mode.  The  second¬ 
mode  lobe-separation  angle  6pp  was  chosen  for  calculations 
for  better  accuracy. 

An  interesting  feature  of  this  laser  is  the  cavity  selection 
of  modes  at  an  intermediate  temperature  of  38  "C,  where  Fig. 
4  suggests  both  lower-order  modes  may  exist  during  a  cur¬ 
rent  pulse.  In  Figs.  8(d)  and  8(e)  the  far  field  is  filtered  with 
the  monochromator  and  the  output  of  the  photomultiplier  is 
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FIG.  3.  Lasing  wavelength  vs  heat-sink  temperature  of  PL174-1-31. 


fed  into  the  sampling  scope  with  a  50-/2  termination.  The  top 
pulse  of  each  picture  is  the  current  pulse  while  the  lower  is 
the  optical  pulse  of  the  laser.  It  was  established  during  the 
experiment,  from  output  pulses  as  in  Fig.  8  and  from  near¬ 
field  pictures,  that  the  laser  switched  modes  during  the  cur¬ 
rent  pulse.  The  discussion  to  follow  explains  this  multimode 
switching  with  current  magnitude  at  an  intermediate  tem¬ 
perature.  Theoretical  results  accompanying  the  experiment 


REFRACTIVE  INDEX  % 

3.600  3598  3596 


FIG.  7.  Experimental  (solid  linel  vs  heat-sink  temperature  and  five-layer 
waveguide  calculated  (dashed  line)  vs  active-region  index. 
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FIG.  8.  Time  display  of  input  laser  current  and  optical  output.  In  (a)  an 
input  current  pulse  is  drawn.  7lh  and  J,hi  are  the  threshold  current  densities 
of  the  first  two  modes,  respectively.  In  (b),  optical  flux  0,  at  the  fundamen¬ 
tal  mode  wavelength  is  shown  to  exist  in  regions  I  and  III.  In  (c),  the  second 
mode  optical  flux  <P2  is  detected  in  region  II  at  the  second  mode  wavelength . 
In  (d),  the  experimental  output  pulse  is  shown  in  the  lower  trace  at  the  fun¬ 
damental  mode  wavelength.  The  top  pulse  is  the  input  current.  In  (e),  the 
lower  trace  is  the  second-mode  output  at  its  wavelength. 


show  the  cavity  selection  is  sensitive  to  gain/mode  profile 
coupling  and  index  change  with  excess  carrier  injection. 


IV.  DISCUSSION 

The  data  from  Fig.  3  were  used  with  the  data  from  Ta¬ 
ble  I  for  a  multilayer  waveguide  analysis  of  the  transverse 
modes.  This  method  was  used  first  to  check  the  accuracy  of 
material  parameters.  In  Fig.  6,  the  dashed  line  shows  a  close 
fit  with  the  solid  measured  curve  of  the  far  field.  The  pip  on 
one  shoulder  is  a  leaky  wave  coupling  to  the  substrate  layer. 
As  previously  mentioned,  Snell’s  law  dictates  a  fixed  refrac¬ 
tive  index  for  the  substrate.  It  was  noticed  in  the  experiment 
that  the  pip  remained  stationary  in  angular  displacement 
with  temperature  change  20-40  °C  and  with  current  drive. 
Also,  the  pip  remains  about  4  dB  below  the  main  lobes  over 
temperature  and  current  variations.  This  means  the  thick¬ 
ness  and  refractive  index  of  layer  3  in  Fig.  2  should  be  fixed  in 
order  to  hold  the  pip  stationary.  The  only  feature  to  change 
for  the  second  mode  is  the  angular  separation  of  the  major 
lobes  designated  as  0pp.  These  data  are  represented  in  Fig.  7 
and  will  be  discussed  subsequently. 

The  only  way  to  match  the  8pp  change  with  temperature 
is  to  slightly  decrease  the  active-layer  refractive  index  with 
increased  current  drive.  The  increase  in  threshold  current  of 
a  mode  with  temperature  occurs  due  to  a  modal  confinement 
change,  internal  quantum-efficiency  decrease,  or  a  carrier- 
lifetime  or  absorption-coefficient  change.20  Near-field  mea¬ 
surements  did  not  show  a  shift  of  the  transverse  mode  with 
temperature.  The  refractive-index  change  of  each  layer  with 
temperature  STj/STcun  be  shown  tobe  1.5X  10  '/°K.:1  For 
a  twenty-degree  range,  20-40  °C,  t577~0.003.  On  the  other 
hand,  the  index  change  due  to  dispersion21  for  a  layer  is  Si]/ 

0.8/eV.  For  SA  =  20  A  from  Fig.  5  (the  first  high-order 
-••  de),  it  can  be  shown  that  Si]^  —  0.003/A. 

The  dispersion  and  heating  effects  on  refractive  index  of 
each  layer  tend  to  cancel.  Therefore,  with  the  evidence  that 
the  near  field  remains  constant  with  temperature,  the  pip  on 
the  shoulder  of  the  far  field  remains  stationary  and  disper¬ 
sion  and  heating  effects  tend  to  cancel,  it  is  assumed  that  the 
dominant  refractive-index  change  is  that  of  the  active  layer 


TABLE  I.  Typical  values  for  ihe  physical  constants  used  in  the  ambipolar 
diffusion  solution .* 


Layer  1 

Layer  2 

Layer  3 

Doping 

p-10  " 

p- 5>  I0r 

n- 5  ■  llV 

A1  concentration 

20% 

0 

2% 

V 

3.48 

variable 

3.5865 

Absorption 

5 

20.8 

2  4 

Width 

2.65  x  10  4 

/*» 

3132 

3649 

3649 

138 

167 

167 

D„ 

81.0 

198 

128 

Dp 

3.67 

4.65 

4.33 

T 

7.69x10  * 

3.59  >  10  4 

8.46  >.  10  4 

D • 

80.6 

15.6 

6.8 

L  * 

7.87x10  4 

2  37x10  4 

2.4  v  10  4 

*  All  quantities  are  in  cgs  units,  with  the  exception  of  mobility  w  hich  is  in  the 
mixed  units  cm/V  sec  and  absorption  which  is  in  cm  '. 
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TABLE  II.  Calculated  values  for  threshold  current  density,  modal  absorp¬ 
tion  and  internal  differential  quantum  efficiency  obtained  from  the  ambtpo- 
lar  diffusion  solution 


Mode  1 

Mode  2 

Threshold  current  density  |kA/cm'l 

15.58 

16.93 

Modal  absorption  icm  '| 

19  74 

6.58 

Internal  differential  quantum  efficiency  (95-) 

60.8 

72.8 

due  to  injection  level.  The  data  of  Tables  I  and  II  and  Fig.  3 
are  used  to  find  the  solution  to  the  ambipolar  diffusion  equa¬ 
tion.  A  rise  in  temperature  is  represented  by  a  reduction  in 
refractive  index  as  required  in  matching  the  far-field  pattern 
of  Fig.  6. 

With  the  relevant  material  parameters  listed  in  Fig.  3,  a 
solution  was  found  for  the  ambipolar  diffusion  equation.  In 
Fig.  9  the  threshold  current  density  for  the  first  two  modes  is 
plotted  as  a  function  of  the  active-region  refractive  index. 
For  larger  rj2  values  (corresponding  to  a  larger  index  step) 
the  second  mode  has  the  lowest  threshold  due  to  superior 
coupling  to  the  excess  carrier  distribution.  But  as  the  active- 
region  index  is  lowered  and  the  second  mode  approaches 
cutoff,  its  coupling  efficiency  diminshes  and  the  fundamen¬ 
tal  mode  emerges  with  the  lowest  threshold. 

Typical  material  parameters  calculated  (fundamental 
mode  at  77-  =  3.5948)  are  shown  in  Table  I.  Since  many  of 
the  parameters  depend  on  injection  levels,  the  quantities  in 
Table  I  vary  slightly  over  the  graphs  presented.  With  these 
calculated  values,  the  usefulness  of  the  ambipolar  equation 
can  be  shown  for  the  active  layer.  The  calculated  hole-injec¬ 
tion  efficiency  is  99  4%  and  the  electron-injection  efficiency 
is  86<?e  for  this  case  This  calculation  predicts  a  crossover 
point  at  a  slightly  lower  r)2  value  than  the  multilayer  mode 
range  in  Fig.  7.  This  is  because  in  the  multilayer  model,  the 
presence  of  the  substrate  with  a  relatively  high  refractive 
index  increases  the  active-region  eigenvalue.  For  instance,  if 
the  buffer-layer  thickness  d ,  were  decreased,  the  active-lay¬ 
er  eigenvalue  would  shift  even  more. 
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FIG.  9.  Ambipolar  diffusion  solution.  Threshold  current  density  for  the 
first-and  second-order  modes  is  plotted  against  the  active-region  index  17- 
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Figure  9  is  useful  in  interpreting  the  phenomena  shown 
in  Fig.  4.  The  fundamental  mode  does  not  occur  below  30  °C. 
At  20  °C,  the  second  mode  has  a  threshold  current  density  of 
about  12  kA/cm'.  However,  abruptly  above  30  °C,  the  fun¬ 
damental  mode  appears.  From  Fig.  8,  a  slight  change  in  tj2 
on  the  order  of  Si 7,  =-  -  0.002  is  sufficient  to  cause  a  cross¬ 
over  in  threshold  currents  between  the  two  modes.  Both 
Figs.  4  and  9  show  that  the  second-mode  threshold  increases 
sharply  after  the  crossover.  Threshold  current  densities  pre¬ 
dicted  by  the  ambipolar  solution  are  in  fair  agreement  with 
experimental  values. 

Modal  output  power  as  a  function  of  current  density  is 
plotted  in  Fig.  10.  Because  the  linear  theory  of  gain  is  as¬ 
sumed,  only  one  mode  at  a  time  is  calculated.  At  power  levels 
above  threshold,  the  modal  gain  is  assumed  pinned  at  its 
threshold  value.  Note  that  the  second  mode,  despite  higher 
threshold  current  density  and  reduced  coupling  efficiency, 
has  a  higher  internal  differential  quantum  efficiency  than  the 
fundamental  mode.  Thus,  for  a  sufficient  overdrive  beyond 
threshold,  the  second  mode  will  eventually  be  the  more  effi¬ 
cient  overall  and  will  dominate.  An  explanation  for  this 
higher  internal  differential  quantum  efficiency  can  be  seen  in 
Table  II.  The  modal  absorption  of  the  second  mode  is  low 
because  a  significant  portion  of  its  energy  is  propagating  out¬ 
side  the  lossy  active  region  whereas  most  of  the  fundamental 
mode  is  well  confined  to  the  active  region.  Since  moda!  loss 
dominates  slope  efficiency,  the  second  mode  has  the  higher 
value.  This  higher  value,  in  conjunction  with  a  higher 
threshold  due  to  poorer  coupling  efficiency,  causes  the  mode 
efficiency  crossover  point. 

With  these  results,  the  phenomena  of  mode  crossover  in 
this  laser  can  be  understood.  At  20  °C  the  second  mode  has 
the  lower  threshold  and  higher  efficiency  and  consequently 
only  the  second  mode  appears.  At  an  intermediate  tempera¬ 
ture.  the  second  mode  will  appear  at  high-enough  current 
overdrive  due  to  better  slope  efficiency.  Above  42  °C,  the 
threshold  for  the  second  mode  becomes  prohibitive.  Figure  8 
depicts  a  time  resolution  between  the  two  modes  at  an  inter¬ 
mediate  temperature.  Figure  8(a)  represents  an  input  drive 
pulse.  The  two  threshold  current  densities  for  the  two  modes 
are  labeled  7,  and  7,  As  the  current  pulse  rises  above 7,  the 
first  mode  turns  on  in  region  I.  Output  flux  <£,  for  the  first 
mode  is  represented  in  Fig.  8(b).  As  the  input  current  pro¬ 
gresses  to  7,,  region  II  is  entered.  The  second  mode  turns  on, 
as  shown  in  Fig.  8(c),  as  <P2.  Finally,  with  the  drop  of  input 
current,  region  III  is  entered  and  first  mode  reappears.  As 
shown  in  Figs.  8!d|  and  8(e),  this  situation  is  observed  experi¬ 
mentally.  In  Fig.  8ld),  the  monochromator  is  set  on  the  first¬ 
mode  wavelength  and  the  lower  pulse  is  the  output  power. 
Reviewing  the  two  output  pulses,  it  can  be  seen  the  two 
modes  have  a  time  separation  during  a  current  pulse.  This  is 
explained  by  the  fact  that,  from  Fig.  9,  the  second-mode 
threshold  can  still  be  reached  slightly  after  crossover.  From 
Fig.  10,  the  second  mode  clearly  is  predicted  to  dominate  the 
first  mode  if  threshold  can  be  reached,  due  to  better  effi¬ 
ciency. 

Theoretical  mode  and  excess  carrier  coupling  for  the 
laser  is  shown  in  Fig.  1 1,  where  the  excess  carrier  profile  is 
plotted  across  the  active  layer  for  the  two  modes  operating  at 
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FIG.  10.  Ambipolar  diffusion  solution. 
Output  power  is  plotted  against  threshold 
current  density  for  the  first  two  modes. 
The  slight  curvature  in  the  mode  1  curve  is 
caused  by  carrier  redistribution.  The  in¬ 
dex  r/2  =  3.5948. 


I 


two  power-output  levels  (threshold  and  0.67  w/tnil.).  Includ¬ 
ed  on  the  carrier  plot  is  the  normalized  modal  distribution. 
Several  factors  should  be  noted.  First,  the  second-mode  dis¬ 
tribution  is  clearly  penetrating  the  passive  n  layer,  indicating 
a  near  cutoff  situation  and  a  modal  loss  dominated  by  the  n- 
layer  absorption.  With  all  curves  normalizd,  a  simple  pro¬ 
duct  of  8p  and  E2,  integrated  over  the  active  region,  shows 
the  fundamental  mode  better  coupled  than  the  second,  as 
expected.  Second,  it  can  be  seen  that  the  carrier  profile  peaks 
on  the  side  where  holes  are  injected.  This  profile  contradicts 
the  notion  that  electron  injection  dominates  and  the  peak 
occurs  on  the  n  side.  Since  the  charge  cloud  is  ambipolar,  it 
must  accommodate  the  low-hole  mobility  to  maintain  quasi¬ 
charge  neutrality.  Hence,  the  profile  is  skewed  toward  the 
hole  injection  or  p  side  of  the  active  layer.  As  expected,  car¬ 
rier  profiles  change  shape  with  current  density. 

As  a  final  note,  we  use  the  refractive-index  dependency 
of  the  active  region  on  temperature  and  hence  current  den¬ 
sity  to  estimate  the  free-ca trier  effect.  The  influence  of  free 
carriers  on  the  refractive  index  in  a  material  is  determined 
from  the  classical  theory  of  dispersion.  The  refractive  index 
in  the  presence  of  free  carriers  tj  is  related  to  the  refractive 
index  in  the  absence  of  carriers  f}  by 


8rjfc  —  K  N,  (17) 

where  AT  is  a  constant  we  wish  to  determine.  By  plotting 
refractive-index  change  versus  carrier  increase,  the  constant 
K  would  be  the  slope  at  a  particular  point 

d(8y(c)  =  Kd{N)  (18) 


or 


*  = 


d(8yfc) 

d{N) 


(19) 


V  m*a)rjJ 


(15) 


Here,  m*  is  the  effective  mass  of  the  carrier,  q  is  the  charge  of 
an  electron,  N  is  the  carrier  density,  eo  —  lire/ A,  and  e0  is  the 
vacuum  permittivity.  Since  the  change  in  index  is  much 
smaller  than  17,  Eq.  ( 1 5)  can  be  rearranged  to 

ItrNq2 


&Vfc  — 


m*eo2y 


(16) 


A  change  in  a  particular  parameter,  say  lifetime  r,  with  injec¬ 
tion  would  be  reflected  in  Eq.  ( 1 7)  having  a  nonzero  second 
derivative.  We  take  8i]k  from  Fig.  7  and  conclude  that  this 
index  change  is  very  nearly  a  linear  function  of  temperature. 
However,  we  find  from  Fig.  4,  carrier  concentration  N,  relat¬ 
ed  to  the  peak  threshold  current,  is  not  quite  a  linear  function 
with  temperature.  Therefore,  under  the  assumption  the  free- 
carrier  effect  is  the  dominant  cause  of  observed  pattern 
change ,  the  constant  K  would  be  precisely 

K_  d  (8i)(Q )/ dT 
d(N)/dT 

This  expression  contains  the  information  that  K  is  not  a  true 
“constant”  but  is  itself  a  slight  function  of  injection  level. 
Nevertheless,  it  is  instructive  to  estimate  K  using  Eq.  (19)  to 
compare  with  conventional  theory.  Using 

qd2 

and  a  straight  line  approximation  to  d  ( N )  over  the  tempera¬ 
ture  range  20-42  °C  we  get 

r  =  3  nsec 


K  =  -  3.1  X  lO”2'  cm3. 


In  other  words,  the  change  in  index  due  to  free  carriers  is 
directly  proportional  to  the  concentration  of  free  carriers 


and 


A=  -  6.2x10 


cm 


t  —2  nsec. 
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FIG.  11.  Excess  carrier  densities  plotted 
across  the  active  region  for  the  ambipo- 
lar  diffusion  equation.  In  |a),  the  second¬ 
mode  distribution  is  plotted.  In  (b),  the 
first-mode  distribution  is  shown.  Note  in 
(a)  that  the  second  mode  “tail"  extends 
into  the  n  cladding. 


These  results  can  be  compared  to  K  =  —  4x  10 —  2 1  cm' 
from  Thompson22  and  K  =  —  (5-ll)X  10~21  cm3  from  Sol¬ 
way.23  Thompson  used  r  =  3  nsec  for  his  calculation.  Sever¬ 
al  other  values  for  K  have  recently  been  reported.  Olsson  and 
Tang24  report  K  =  —  (4.9  ±  0.4)x  10-21  for  r  =  2.2  nsec. 
Henry  et  a/.1*  report  K  -  —  (1.8  ±  0.4)  X  10  21  for  lasers 
with  carrier  densities  estimated  at  1.02  X  10IRcm~\ 

V.  CONCLUSION 

A  solution  for  the  ambipolar  diffusion  was  presented 
for  the  active  region  of  a  double  heterojunction  laser.  Along 
with  the  solution  of  modes  from  Maxwell's  equations  in  a 
multilayer  waveguide,  this  combined  description  was  ap¬ 
plied  to  a  broad-area  heterojunction  laser.  The  refractive- 


index  change  of  the  active  region  with  temperature  was  esti¬ 
mated  from  far-field  pattern  matching  of  the  transverse 
modes.  Transverse-mode-selection  conditions  were  predict¬ 
ed  with  the  diffusion  equation  and  compared  with  the  trans¬ 
verse-mode  selection  of  the  sample  laser. 

Mode  threshold  was  shown  to  depend  on  coupling  effi¬ 
ciency  between  excess  carriers  and  modal  distribution  in  the 
active  layer.  The  differential  internal  quantum  efficiency  de¬ 
pends  primarily  on  mode  absorption  and  injection  effi¬ 
ciency.  Mode  crossover  occurs  when  the  mode  with  the 
highest  threshold  also  has  the  highest  slope  efficiency. 

Carrier  transport  in  a  semiconductor  laser  is  an  ambi¬ 
polar  diffusion  process  which  can  be  described  using  the  re¬ 
laxation  time  approximation.  In  bulk  regions,  quasi-charge 
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neutrality  is  satisfied  because  of  the  strong  internal  field 
between  electrons  and  holes.  The  active  region  of  a  semicon¬ 
ductor  laser  undergoes  double  injection:  hole  injection  from 
the  p  side  and  electron  injection  from  the  n  side.  In  devices 
with  good  injection  efficiencies,  the  injected  carrier  profile  is 
higher  at  the p-p,  rather  than  at  the p-n.  interface  because  of 
low  hole  mobility. 

The  device  modeled  here  shows  a  distinct  peak  of  carri¬ 
ers  on  the  p-p  interface.  The  slope  of  the  carrier  profile  in¬ 
creases  with  increased  current  density.  Since  these  carrier 
profiles  demonstrate  significant  skewing,  it  is  clear  that  this 
can  affect  mode  selection  if  the  modal  shape  couples  poorly 
with  the  carrier  concentration. 
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